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Abstract. 

The space of parameters characterizing an inflationary primordial power spectrum 
with small superimposed oscillations is explored using Monte Carlo methods. The 
most interesting region corresponding to high frequency oscillations is included in 
the analysis. The oscillations originate from some new physics taking place at the 
beginning of the inflationary phase and characterized by the new energy scale M c . It 
is found that the standard slow-roll model remains the most probable one given the 
first year WMAP data. At the same time, the oscillatory models better fit the data 
on average, which is consistent with previous works on the subject. This is typical of 
a situation where volume effects in the parameter space play a significant role. Then, 
we find the amplitude of the oscillations to be less than 22% of the mean amplitude 
and the new scale M c to be such that H/M c < 6.6 x 1CP 4 at la level, where H is the 
scale of inflation. 



PACS numbers: 98.80.Cq, 98.70.Vc 
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1. Introduction 

The possibility that the Cosmic Microwave Background (CMB) data contain unexpected 
features which would be signatures of non standard physics has recently been studied 
in Refs. [H E] . These features generally consist of superimposed oscillations in the CMB 
multipole moments Ce. They could originate either from trans-Planckian effects taking 
place in the early universe El HI [HI El [3 El O [101 CU [121 [EH [HI [13 HE1 HZl UHl HH1 1201 
EH E21 ESI EH ESI EE1 E3 EHl EH1 EDI EH E2| or from other phenomena, for instance of 
the type of those investigated in Refs. [331 E3 E3 EHl 133 EE] ■ 

In Refs. jU E] , it has been demonstrated that the presence of superimposed 
oscillations can improve the fit to the first year Wilkinson Microwave Anisotropies Probe 
(WMAP) data [3J1 HOI EJ H2 EH , see also Refs. jmiSj. The corresponding drop in the 
X 2 was found to be A% 2 ~ 10. Then, based on the so-called F-test, it has been argued 
that this drop is statistically significant. However, it should be clear that the study 
carried out in Refs. Ej has only provided a hint with regards to the possible presence 
of non trivial features in the multipole moments. In order to go further, it is necessary 
to explore the parameter space in details, taking into account, of course, the parameters 
characterizing the oscillations, namely the amplitude, the frequency and the phase. This 
task is delicate because the numerical computation of the multipole moments in presence 
of rapid oscillations in the initial power spectrum necessitates an accurate estimation 
of the CMB transfer functions and of the line of sight integrals which, in turn, requires 
to boost the accuracy with which the CVs are computed (here by means of the CAMB 
code jin])- This causes an important increase of the computational time and instead 
of a few seconds, the computation of one model now takes a few minutes. Therefore, 
the exploration of the full parameter space (i.e. the one including the cosmo logical 
parameters characterizing the transfer functions and the primordial parameters encoding 
the shape of the initial power spectra) remains beyond the present day computational 
facilities capabilities, at least in a reasonable time. For this reason, in this paper, we 
restrict ourselves to the exploration of the so-called fast parameter space only, i.e. the 
space of parameters characterizing the primordial spectrum, including of course the 
oscillatory parameters. These parameters are usually called "fast" because restricting 
the exploration to the corresponding parameter subspace requires the calculation of the 
transfer functions only once, hence a huge gain in term of computational time. However, 
in the present work, the correct computation of the superimposed oscillations that are 
transfered from the power spectrum to the multipole moments requires a very accurate 
evaluation of the line of sight integrals. Therefore, even in the fast parameter subspace, 
each determination of the CVs remains a heavy time consuming task PIE]- The main 
new result of the present article is that, for the first time, the shape of the likelihood 
function in this space including, and this is a crucial point, the region containing the 
best fit found in Refs. [TJ E] is presented. 

The present study also allows us to gain some new insight into the statistical 
significance of the superimposed oscillations. In particular, we derive the one- 
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dimensional marginalized probabilities and the mean likelihoods for the oscillatory 
parameters. Our result indicates that the slow-roll models are still the most probable 
ones. However, the marginalized probability distribution of the amplitude possesses 
a long tail originating from the fact that the mean likelihood is peaked at values 
corresponding to models with oscillations. From the above-mentioned distributions, 
we derive constraints on the parameters characterizing the oscillations. 



Finally, some comments are in order about the previous works on the subject |191 

EElEBj- In those references, the shape of the likelihood function in the parameter space, 
including cosmological parameters, was studied and its hedgehog shape highlighted. 
However, the exploration was limited to the low frequency models in order to avoid 
the above-mentioned numerical difficulties. As a consequence, the region where the fit 
found in Refs. [H |2] lies was not considered. As mentioned above, the exploration of 
this region is a crucial ingredient for the present article and plays a role in constraining 
the oscillatory parameters. 

This paper is organized as follows. In the next section we briefly recall the 
theoretical model used to generate primordial oscillations, together with the numerical 
implementations. In Sec. El marginalized probabilities and constraints on the fast 
parameters are discussed. We give our conclusions in Sec. |U 

2. The fast parameter space 

We assume that inflation of the spatially flat Friedman-Lemaitre-Robertson- Walker 
(FLRW) space-time is driven by a single scalar field (f>(t). Then, the scalar fluctuations 
can be characterized by the curvature perturbations ( and the tensor perturbations by 
a transverse and traceless tensor hff ^7J EH1 HH EH]- If the trans-Planckian effects 
are taken into account, the corresponding power spectra take the form of a standard 
inflationary part plus some oscillatory corrections. The standard part is parameterized 
by the Hubble parameter during inflation H, the slow-roll parameters e 1 and e 2 , and the 
pivot scale k* [HJ HHl EH EH] • The two slow-roll parameters that we use are defined by 
e 1 = —H/H 2 and €2 = 2(e 1 — 5) with 5 = —<p/{H(j)). The oscillatory part is determined 
by three new parameters. Two of them are the modulus \x\ and the phase <p of a complex 
number x which characterizes the initial conditions and the other is o~ = H/M c , M c 
being a new energy scale. Note that o~ is evaluated at a time T] during inflation which 
is a priori arbitrary but, and this is the important point, does not depend on k. In the 
following we will choose this time such that k*/a = M c where a is the scale factor at 
time 7] . Explicitly, for density perturbations, the power spectrum reads PQI2] 



and a similar expression for the gravitational waves. Typically, the amplitude of the 
oscillations is given by \x\a Q while the wavelength is such that Ak/k = 7ra /e l . As 
discussed in Ref . pQ , the factor x parameterizing the initial conditions plays an important 
role in our analysis. The quantum state in which the cosmological fluctuations are 




k 3 P ( = k 3 P s J 1 - 2|a;|a cos & 



In(£)+V] + O(\x\a e 2 ) 



(1) 
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"created", when their wavelength equals the new characteristic scale M c , is a priori 
not known. Usually, it is assumed that this quantum state is characterized by the 
choice x — 1. In the present article, following Refs. [2HJ 1221 EH Elj , we generalize these 
considerations and describe our ignorance of the initial state by the parameter x which 
is therefore considered as a free parameter. From the point of view of the statistical 
study presented here, this turns out to be very important. Indeed, as mentioned in 
Ref. [T] [see the discussion after Eq. (8)], this allows us to decouple the amplitude of the 
oscillations from the value of the frequency which, therefore, can now be considered as 
independent parameters. As a consequence, high frequency waves no longer necessarily 
have a small amplitude as it is the case if one chooses x — 1. Let us also notice that, 
despite the fact that the parameterization of the initial conditions used is more general, 
we have nevertheless restrict our considerations to the case where x is scale independent 
over the scales of interest. 

In order to compute the resulting CMB anisotropies, we use a modified version of 
the CAMB code [46J whose characteristics and settings are detailed in Refs. [UEj- The 
exploration of the parameter space is performed by using the Markov Chain Monte- 
Carlo methods implemented in the COSMOMC code [HI] together with our modified CAMB 
version. As already mentioned, in order to avoid prohibitive computational time, only 
the "fast parameter space" is explored. Notice again that, even in this case, the required 
accuracy (see Refs. PJE]) renders the numerical computation much longer than the one 
necessary to the exploration of the "full parameter space" in the standard inflationary 
case. 

Since we restrict ourselves to the "fast parameter space", we have to fix several 
degrees of freedom. Firstly, the cosmological parameters determining the CMB transfer 
functions should be chosen. Since the oscillatory features in the multipole moments 
are small corrections to the standard inflationary predictions, a reasonable choice is 
to fix the ("non primordial") cosmological parameters to their best fit values obtained 
from the vanilla slow- roll power spectrum, i.e. from a pure slow-roll model without 
any additional feature, in a flat ACDM universe PHEJE2]- Indeed, it has been shown 
in Ref. |2E3 that the oscillatory parameters and cosmological parameters are not very 
degenerate. This amounts to take h = 0.734, Vt h h 2 = 0.024, tt dm h 2 = 0.116, fi A = 0.740, 
r = 0.129, z Te = 14.3 and lO 2 ^ = 1.049, 9 being approximately the ratio of the sound 
horizon to the angular diameter distance, see Ref. [5T] . 

Secondly, priors on the remaining fast parameters have to be imposed. For the 
standard inflationary parameters, e 1 , e 2 and the logarithm of the initial amplitude of 
the scalar power spectrum at the pivot scale, lnP sca i ar , wide flat priors around their 
standard values have been chosen. Let us now discuss the three oscillatory parameters. 
For the phase, the sampling has been performed on the quantity ip assuming a flat prior 
in the interval [0, 27r] because there is also a sine function in the term that we have not 
written explicitly in Eq. (JTJ). For the amplitude of the oscillations, i.e. \x\a Q , we have 
also imposed a flat prior and have required this parameter to vary in the range [0, 0.45], 
the upper limit being chosen in order to avoid negative primordial power spectra. 
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Figure 1. ID marginalized probabilities (solid curve) and normalized mean likelihoods 
(dashed curve) for the slow-roll parameters e ± , e 2 and the scalar amplitude P S caiar, m 
the case of the reference model. The three-dimensional plots show the normalized 
mean likelihoods and the corresponding ler and 2a contours of the 2D marginalized 
probabilities. 

The case of the last oscillatory parameter, i.e. the frequency, should be treated 
with great care. The frequency of the superimposed oscillations is given by 2e 1 / a a and, 
hence, we have chosen to sample the quantity log (l/a ). At this point, it is necessary to 
clarify what we exactly mean by "superimposed oscillations" . Indeed, if the frequency is 
very low and is such that less than half of a period covers the whole observable range of 
multipoles, then the resulting effect is just a modification of the overall amplitude of the 
power spectra and we don't really have oscillations in the multipole moments anymore. 
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In this regime, there is a strong degeneracy between the amplitude and P sca iar- Clearly, 
it is not desirable to reach this extreme case. Another regime corresponds to the case 
where the frequency of the superimposed oscillations is of the same order of magnitude 
than the frequency of the acoustic oscillations. In this situation, new peaks would appear 
in the C/s and this could be problematic since the WMAP data are well described by 
the standard inflationary model. As a matter of fact, the amplitude of this type of 
oscillations is strongly constrained, see e.g. Ref. [21] • Finally, there is the case where 
high frequency oscillations are present. This is the case we are mostly interested in since 
it contains the improved fits found in Refs. PQI2]- Nevertheless, very high values of the 
frequency lead to oscillations from multipole to multipole and any dependency on the 
frequency is in fact lost. It is clear that we would also like to avoid this extreme regime. 
All the above considerations has led us to restrict ourselves to log (l/cr ) G [2.3, 3.8] with 
a flat prior. This is equivalent to imposing a Jeffreys-like prior on o~ . 

Our next step has been to generate Markov chains with COSMOMC for a standard 
slow-roll inflationary model without oscillation (i.e. for our reference, or vanilla, model) 
as well as for an oscillatory model characterized by ip, \x\a and log(l/c ) with the above 
priors, both cases being investigated with the non-primordial cosmological parameters 
chosen before. The standard inflationary model can be used as a reference for the 
restricted framework of the "fast parameter space". For the two models, the Markov 
chains have been started from widespread point in the initial parameter space and 
stopped when the posterior distributions no longer evolve significantly which corresponds 
to about 250000 elements. The generalized Gelman and Rubin i?-statistics implemented 
in COSMOMC j^H EHj ^ s found to be less than 5% for the slow-roll parameters and equal 
to approximately 9% for the oscillatory parameters. However, the phase ip remains 
unconstrained. Let us recall that the ill-statistics gives, for each parameter, the ratio 
of the variance of the chain means to the mean of the chain variances. In other words, 
it quantifies the errors in the parameter distributions obtained from the Markov chains 
exploration. 

3. Results 

In Fig. Q we have plotted the one dimensional marginalized probabilities and the 
normalized mean likelihoods for the slow-roll parameters and the scalar amplitude 
in the case of the reference model. Three-dimensional plots of the mean likelihoods, 
obtained by averaging with respect to one parameter, are also represented with the 
la and 2a contours of the two-dimensional marginalized probabilities which appear 
through the surface. These plots are obviously consistent with the previously derived 
constraints [32], except that fixing the cosmological parameters has made the likelihood 
isocontours slightly tighter. 

In Fig. El we have plotted the ID marginalized probabilities and the normalized 
mean likelihoods for the parameters of the oscillatory model. In Figs. El and EJ the la and 
2a contours of the 2D marginalized probabilities and the normalized mean likelihoods 
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are displayed. They have been plotted in various planes of parameters, the parameters 
pairs being chosen to be the most correlated ones. 

Let us now discuss the constraints that can be put on the amplitude of the 
oscillations, \x\a . As can be seen in Fig. |2J the vanilla slow-roll model, corresponding 
to \x\a = 0, remains the most probable one with the currently available WMAP data. 
From the probability distribution, one finds that the 2a marginalized upper limit is 
given by 

\x\a < 0.11 . (2) 

However, the most striking feature of Fig. 121 is that the corresponding mean likelihood 
does not behave similarly to the marginalized probability, as it is the case for the other 
parameters. Indeed, it exhibits a maximum around \x\a ~ 0.1 which is consistent with 
the best fits found in Refs. QjJ E] . This unusual behavior is expected if volume effects 
are present (more precisely, since we are dealing with quantities defined as integrals over 
the parameter space, the volume effects are important when the corresponding kernels 
are spread over this space with large values). This highlights the existence of strong 
non-linear correlations between \x\a and the other parameters. 

As can be seen in Fig. strong correlations appear between the amplitude and 
the frequency of the oscillations. On the three-dimensional plot [|x|cr , log (l/cr )], one 
remarks that the mean likelihood is non- vanishing either for large values of log(l/<r ), 
when \x\a ^ 0, or in a "long thin tube" centered around |x|cx = 0, the amplitude 
of which being much smaller (recall that A% 2 ~ 10 between these two models). This 
behavior is expected since any frequency is possible provided its amplitude is very small 
hence the "long thin tube" around \x\a = 0. On the other hand, for large values of 
log (1/<t ), one recovers a peak which corresponds to the oscillations fitting the outliers 
at relatively small scales [HI2]- The marginalized probability is peaked at |x|a = 
because the tube has a bigger statistical weight, i.e. occupies a larger volume in the 
parameter space. Note also that this is the very existence of the likelihood peak at large 
log(l/<r ) which explains the long tail of the ID marginalized probabilities for \x\a . 

On the three-dimensional plot [e i; log (1/<t )], the wavelets in the mean likelihood 
originates from the degeneracy existing between a and e 1 , for a fixed value of the 
frequency 2e l /a a . Furthermore, the correlation between e 1 and e 2 52] is still present 
and, as a consequence, a correlation between e 2 and log(l/a ) also exists. We also 
remark that the correlations appearing when the oscillatory parameters are taken into 
account in the analysis seem to slightly tighten the confidence contours on the standard 
inflationary parameters and to cause the appearance of some multimodal features. 

Another interesting property concerns the parameter log (1/<t ). As can be seen in 
Figs. |2l El an d EJ both the mean likelihood and the marginalized probability are peaked 
at a high value of log (l/cr ). Quantitatively, a lcr constraint on the energy scale M c at 
which the new physics shows up can be derived 

ff ° = w c < 6 - 6 x 10-4 • ( 3 ) 

Let us stress that this constraint is valid only if the amplitude and the frequency of the 
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Figure 2. ID marginalized probabilities (solid curves) and normalized mean 
likelihoods (dashed curves) for the oscillatory model. 

superimposed oscillations are considered as independent parameters. Indeed, it mainly 
comes from the region where the likelihood is strongly peaked and, as already noticed, 
this region does not exist if the parameter x is chosen such that x — 1. In this case, the 
likelihood is rather flat and, as a consequence, one would obtain a weaker contraint on 
the ratio H/M c . 

Finally, we have found that the phase i/j is not constrained as is clear from Fig. El 
To end this section, one can ask how well the oscillatory model fits the data on 
average, compared to the standard slow-roll one. This can be estimated by comparing 
the mean likelihood over all the parameters in each of these models [HI] . For the reference 



Exploring the Superimposed Oscillations Parameter Space 



9 




Figure 3. la and 2a contours of the 2D marginalized probabilities and normalized 
mean likelihoods (see also Fig. . 

model one gets ln£ sr ~ —715.1 while the oscillatory model has ln£ wig ~ —712.9, 
hence a ratio of exp (715.1 — 712.9) ~ 9 in favor of the oscillatory model. This shows 
that including the oscillatory parameters permits to improve the goodness of fit on 
average [ST] . 
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Figure 4. la and 2a contours of the 2D marginalized probabilities and normalized 
mean likelihoods (see also Fig. ■ Strong correlations appear between the amplitude 
and the frequency of the oscillations (on the left). 

4. Discussion and Conclusion 

In this paper, we have explored, by means of Monte Carlo methods, the fast 
parameter space of an inflationary model with superimposed oscillations. This restricted 
framework is, at the time of writing, the only way to derive constraints in a reasonable 
computational time if the region where high frequency oscillations that better fit the 
CMB anisotropics outliers at relatively small scales [H E| is included. Among the 
main results derived in the present article are two constraints, one on the amplitude, 
\x\a < 0.11, and the other of the energy scale M c , H/M c < 6.6 x 10~ 4 . 

The overall situation is quite interesting: on one hand, the fact that the 
marginalized probability is peaked at a value corresponding to a vanishing amplitude 
shows that the most probable model remains the standard slow-roll one for which 
\x\a = 0. On the other hand, this distribution exhibits a long tail in the regions 
corresponding to non-vanishing amplitudes of the oscillations, precisely where the mean 
likelihood function is peaked, i.e. around \x\a ~ 0.1. Moreover, the ratio of the vanilla 
slow-roll model to the oscillatory model total mean likelihoods is about 9, in favor of 
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the oscillatory model. This shows that, on average, the oscillatory model fits better the 
first year WMAP data. 

The interpretation of the situation described above is as follows. The marginalized 
probabilities are quantities which are especially sensitive to "volume effects", i.e. to the 
shape of the likelihood in regions of the parameter space where this one is significant. On 
the other hand, the mean likelihood is sensitive to the absolute value of the likelihood 
regardless to its occupied "volume". Therefore, if in the parameter space there are 
regions which, at the same time, correspond to good fits and occupy a quite confined 
volume, then these regions will appear particularly significant from the mean likelihood 
point of view but will be "diluted" and, hence, will appear less significant from the 
marginalized probabilities point of view. Ultimately, the marginalized probability is 
the relevant quantity, that is to say is the quantity which should be used in order to 
conclude about the statistical meaning of a given model [Slj. It will be interesting to 
study how the situation evolves when better data become available, in particular it will 
be interesting to see if the models with oscillations can occupy a volume which makes 
them probable from the marginalized probabilities point of view. This is one of the 
reasons why the next WMAP data release will be of great interest. 
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